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In this paper we show the series of Greenberger-Horne-Zeilinger paradoxes for N maximally 
entangled iV-dimensional quantum systems. 
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I. INTRODUCTION 

The Greenberger-Horne-Zeilinger (GHZ) correlations, discovered in 1989, started a new chapter in the research 

related to entanglement. To a great extent this discovery was responsible for the sudden renewal in the interest 

| in this field, both in theory and experiment. All these developments finally led to the first actual observation of 

three qubit GHZ correlations in 1999 [EJ, and as a by product, since the experimental techniques involved were 

of the same kind, to the famous teleportation experiment [Q]. 

In j^, |j it has been shown that GHZ paradox can be extended to the correlations observed in quantum 

systems consisting of N + 1 maximally entangled TV-dimensional quantum objects (so called quiVits), where 

./V is an arbitrarily high integer number. The existence of the GHZ paradox for three maximally entangled 

three dimensional quantum systems has been shown in . It is worth mentioning that two different approaches 

have been presented in to derive GHZ paradoxes. In || the series of paradoxes has been derived for 

the correlation function that results from the correlations observed in so called unbiased symmetric multiport 

■ beamsplitters (for the description of such devices see below) whereas Q presents derivation based on the relations 

| between operators. 

In this paper we would like to show that GHZ-type paradoxes exist also in the case of correlations expected 

in gedanken experiments involving N maximally entangled quiVits. This is an important question as the 

existence of such paradoxes allows us to deeper understand the structure of quantum entanglement, which is the 

fundamental resource in quantum information as well as the nature of non-classicality of quantum correlations 

(by non-classicality we understand lack of local realistic description of such correlations; the alternative measure 

fi i| of non-classicality can be a measure based on the notion of non-separability). 

To this end, we shall study a GHZ-Bell type experiment in which one has a source emitting N qu/Vits in a 

specific entangled state of the property, that the quTVits propagate towards one of N spatially separated non 

conventional measuring devices operated by independent observers. Each of the devices consists of an unbiased 

qh symmetric multiport beam splitter || (with N input and N exit ports), N phase shifters operated by the 

observers (one in front of each input), and N detectors (one behind each exit port). 
»" ■ 

•jj ■ II. UNBIASED MULTIPORT BEAMSPLITTERS 

a, 

An unbiased symmetric 2iV-port beam splitter is defined as an TV-input and TV-output interferometric device 
which has the property that a beam of light entering via single port is evenly split between all output ports. 
I.e., the unitary matrix defining such a device has the property that the modulus of all its elements equals -^=. 

An extended introduction to the physics and theory of such devices is given in [[|. Multiport beam splitters 
were introduced into the literature on the EPR paradox in |Io|-|l2| in order to extend two qubit Bell-phenomena 
to observablcs described as operators in Hilbcrt spaces of the dimension higher than two. In contradistinction 
to the higher than 1/2 spin generalizations of the Bell-phenomena [fl3|-[l9||, this type of experimental devices 
generalize the idea of beam-entanglement [^o] . Unbiased symmetric multiport beam splitters are performing 
unitary transformations between "mutually unbiased" bases in the Hilbert space p3-|2(|. They were tested in 
several recent experiments and also various aspects of such devices were analyzed theoretically p9| , ^o| . 

We shall use here only multiport beam splitters which have the property that the elements of the unitary 
transformation which describes their action are given by 

TjN _ J_ ( m -l)(m'-l) m 

u m,m' — /J^^N ' V ' 



where jjy = exp(i^) and the indices m, ml denote the input and exit ports. Such devices were called in 
the Bell multiports. 
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III. QUANTUM MECHANICAL PREDICTIONS 



Although in this paper we consider only the situation in which there is an equal number of observers and 
input ports in each multiport beamsplitter it is instructive to derive necessary formulas for the more general 
case in which there are M observers each operating 2N port beamsplitter. The initial M quNit state has the 
following form: 

N M 

iv>(M)}=-=^ni m ><> ( 2 ) 

v m— 1 1=1 

where \m)i describes the l-ih qaNit being in the m-th beam, which leads to the m-th input of the Z-th multiport. 
Please note, that only one quTVit enters each multiport. However, each of the quiVits itself is in a mixed state 
(with equal weights), which gives it equal probability to enter the local multiport via any of the input ports. 

The state (|j) seems to be the most straightforward generalization of the GHZ states to the new type of 
observables. In the original GHZ states the number of their components (i.e., two) is equal to the dimension of 
the Hilbert space describing the relevant (dichotomic) degrees of freedom of each of the quiVits. This property 
is shared with the EPR-type states proposed in || for a two- multiport Bell- type experiment. In this case the 
number of components equals the number of input ports of each of the multiport beam splitters. We shall not 
discuss here the possible methods to generate such states. However, we briefly mention that the recently tested 
entanglement swapping ]3l]-|33|] technique could be used for this purpose. 

As it was mentioned earlier, in front of every input of each multiport beam splitter one has a tunable phase 
shifter. The initial state is transformed by the phase shifters into 

JV M 

mi) 1 ) = -/=£ HeM^DHi, (3) 

v m=l 1=1 

where 0™ stands for the setting of the phase shifter in front of the m-th port of the Z-th multiport. 

The quantum prediction for probability to register the first photon in the output k\ of an 2N - port device, 
the second one in the output k% of the second such device ,. . ., and the M-th one in the output Um of the M-th 
device is given by: 

PQ\i{ki, . . . , k M |0i, ■ • ■ , 4>m) = 
(^) J,/+1 IE™= 1 exp(zE J = 1 C)nf= 1 7^ 1)(fc "- 1) | 2 = 



= [N + 2J2Z >m ,co S (zZi*KC')} ■ ( 4 ) 



where A^J 7 ^ = 0™ — <fi" 1 + — l)(m — to'). The shorthand symbol 4>k stands for the full set of phase 

settings in front of the k-th multiport, i.e. (fy, <jj^, , . . , (f>% . 

To efficiently describe the local detection events let us employ a specific value assignment method (called 
Bell number assignment; for a detailed explanation see again &), which ascribes to the detection event behind 
the m-th output of a multiport the value 7^ _1 , where 7jv = exp(i^). With such a value assignment to the 
detection events, the Bell-type correlation function, which is the average of the product of the expected results, 
is defined as 

E(<f>x, ■ ■ -,<Pm) = 

= EL-Ms=i IL=i lN~ lp (ki, • • • , k M \&,- ■ ■ , fu) (5) 

and as we shall see for the quantum case it acquires particularly simple and universal form (which is the main 
purpose for using this non-conventional value assignment). 

The easiest way to compute the correlation function for the quantum prediction employs the mid formula of 
(§: 

Eqm(4>i,- ■ - ,4>m) = 

ttM (fci-l)(m-m'+l) _ 

x Il;=i In — 
- (i?) M+1 ZZ, M exP (< ELi(C " C')) 

x Il^iEfe=i7^ ! " 1)(m " m ' +1) - (6) 
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Now, one notices that Ylk=i 7w' ™ + ^ differs from zero (and equals to N) only if m — ml + 1 = modulo 

N. Therefore we can finally write: 

Eqm(<Pi, ■ ■ -,<Pm) 

= lvEm = i e WEl=l^ ' )l ( 7 ) 

where <p™' m+1 = 0™ — and the above sum is modulo N, i.e., 0f +1 = 

One can notice here a striking simplicity and symmetry of this quantum correlation function (R). ft is valid 
for all possible values of M (number of quiVits) and for all possible values of N > 2 (number of ports). For 
N = 2, it reduces itself to the usual two qubit, and for N = 2, M > 2 the standard GHZ type multi-qubit 
correlation function for beam-entanglement experiments, namely cos 1 0; 1 ' 2 ) 0- The Bell-EPR phenomena 
discussed in j|] are described by (Q) for M = 2, N > 3. 

Even for N = 2, M — 1 the function @ describes the following process. Assume that a traditional four-port 
50-50 beam splitter, is fed a single photon input in a state in which is an equal superposition of being in each of 
the two input ports. The value of (Q) is the average of expected photo counts behind the exit ports (provided 
the click at one of the detectors is described as +1 and at the other one as —1), and of course it depends on 
the relative phase shifts in front of the beam splitter. In other words, this situation describes a Mach-Zehnder 
interferometer with a single photon input at a chosen input port. For AT = 3, M = 1 the same interpretation 
applies to the case of a generalised three input, three output Mach-Zehnder interferometer described in B, 
provided one ascribes to firings of the three detectors respectively 73 = a = exp(i^), a 2 and a 3 . 

The described set of gedanken experiments is rich in EPR-GHZ correlations (for M > 2). To reveal the 
above, let us first analyze the conditions (i.e. settings) for such correlations. As the correlation function (0) is 
an average of complex numbers of unit modulus, one has \Eqm{4>i, • ■ ■ , 4>m)\ < 1- The equality signals a perfect 
EPR-GHZ correlation. It is easy to notice that this may happen only if 

A I M M 

1=1 1=1 1 = 1 

where k is an arbitrary natural number. Under this condition E(4>i, • • • , 4>m) — This means that only those 
sets of M spatially separated detectors may fire, which are ascribed such Bell numbers having the property that 
their product is 7^. Knowing, which detectors fired in the set of M — 1 observation stations, one can predict 
with certainty which detector would fire at the sole observation station not in the set. 



IV. GHZ PARADOXES FOR N MAXIMALLY ENTANGLED QLWITS. 

Now we show the derivation of GHZ paradoxes for N maximally entangled quiVits. However, first it is 
instructive to consider the simplest case, which is three qutrits. The correlation function for such an experiment 
reads 

Eqm(4>1,4>2,4>3) 

k=l \ 1=1 / 



The experimenters perform three distinctive experiments. In the first run of the experiment we allow the 
observers to choose the following settings of the measuring apparatus, <j>i — <j)2 — (0, j, -^-) — (f>, 03 = (0, 0, 0) = 
</>' in the second run they choose (j>2 = 4>3 — 4>i <f>i = 4>' whereas in the third run they fix the local settings of 
their fritters on <f>i = <f>3 = 4>, 4>2 = 4>' ■ 

Now, let us calculate the numerical values of the correlation function for each experimental situation. Wc 
easily find that for all three experiments this value, due to the special form of the correlation function, is the 
same and reads 

EQM{(j>, 0, <t>') = E QM {<\> ', <t>, (f) = E, 



qm{<P,<P ,<P) 
2tt o 

exp(-t — ) = a , (9) 



i.e., we observe perfect correlations. 
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The (deterministic) local hidden variable correlation function for this type of experiment must have the 
following structure j|] : 

r 3 

EBv(<h,<f>2,<h)= / l[lk(0k,X)p(\)d\. (10) 

jA k=i 

The hidden variable function I k ((f> k , A), which determines the firing of the detectors behind the k-th multiport, 
depends only upon the local set of phases, and takes one of the three possible values a, a 2 , a 3 = 1 (these values 
indicate which of the detectors is to fire), and p(X) is the distribution of hidden variables. The local realistic 
description of the experiment is only possible if the local hidden variable correlation function defined above 
equals that of quantum ones for the given set of phase shifts. 

Because in each experiment described above we observe perfect correlations, i.e., the quantum mechanical 
correlation function takes value a 2 , the local realistic description is possible if and only if one has Ejivifi, §■> < t>') — 
Ehv(<P\ 0, <ft) — E>hv(4>i c / ) 'j ( I ) ) — a2 ■ Taking into account the structure of local hidden variable correlation 
function one obtains the set of three equations 

I 1 ($,\)I 2 &\)I 3 $,\) = a 2 
h($'^)h($A)h{lX)=a 2 

h(l\)I 2 ($',\)I 3 ($ 1 X)=a 2 (11) 
After multiplication of the above equations one arrives at: 

3 3 

l[l k ($',X)l[l k (lX) 2 = (a 2 ) 3 = l, (12) 

fe=l k=l 

which can be also written in the following form 

3 3 

Y[l k ($',\) = l[l k (l\). (13) 

fc=i fe=i 

We have used the property of the hidden variable functions: I k [4>, A) 2 = I k (<fi, A)* for every A. However, because 
Eqm{^\ <t>', </>') = 1 we must also have 



nw'A)=i, (w) 

which, because of (|l3|), gives 

3 

Y[h($,X) = l (15) 

k=l 

for every A. This in turn implies that 

E QM ($JJ) = 1, (16) 

which means that local hidden variables predict perfect correlations for the experiment when all observers set 
their local settings at (f>. However, the quantum prediction is that 

E QM ($JJ) = ~. (17) 

Therefore, we have the contradiction: 1 = — 4. 

This contradiction is of the different type than the one derived in || although it has been obtained in the 
similar way, i.e., the perfect correlations have been used to derive it (equations ( pi] ) and ([l2|)). Here local hidden 
variables imply a certain perfect correlation, which is not predicted by quantum mechanics, whereas in j^] as 
well as in Q the contradiction is that both theories predict perfect correlations but of the different type. 

Now, we employ the above procedure for the case when we have an arbitrary odd number of multiports and 
qu-/Vits N = 2m + 1. As we have seen in in the above derivation the crucial point is to find the proper phases 
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for the multiports. Let us choose for the first gedanken experiment the following ones <f>\ = 02 = 03 = ■■■ = 

(0, 2mTT> 3mTT> • • • ' = I = (0, 0, 0) = 0'. 



02m 



As before, in the next run of the experiment we choose the same phases but we change the role of observers 
such that in the second run the first one chooses 0' while the rest of them choose 0, in the third run the third one 
chooses 0' while the rest of them choose 0, etc. Again, the value of the correlation function for each experiment 
is the same 



Eqm(<P, ■•■,<?') = Eq M (</>', 



1 



2m + 1 



2m exp 



. 2m7r 
'2m + 1 



+ exp i 



E QM (<p, ...,</>', 
4m 2 7r 

' ' = 7at" 



2m + 1 



(18) 



Using (18) and the structure of the hidden variables correlation function, i.e., Ehv{4>1i ■ ■ ■ > 0iv) = 
Ia rifc=i h($k, X)p(\)dX, we arrive at the set of N equations h{$' , X) h($, X) = j N m (i = 1,2,..., N). 

After multiplying them by each other we obtain 



2m+l 



2m+l 



fc=l 



2m+l 

n 

fc=l k=l 
[2m+l 



n h($\x)= n h($,x)- 2m = } [ i k {ix) 



(19) 



for every A. Because Eqm(4>', ■ ■ ■ , = 1 one must have rife"i h{4>\ X) = 1, which according to ( pL9[ ) 
gives IIfc=i" 1 4(0, A) = 1 for every A. Thus, local hidden variables imply the following perfect correlation 
Eqm(4>> •■■i0) = 1? whereas quantum mechanics gives 



E, 



1 



QM( 



2m +1 



2m exp 



.2m+ 1 

-i>x rTi" 

2m + 1 



'.(2m+ l)2m 

-exp | i 7r 

1 [ 2m + 1 



-2m + 1 
2m + 1 



Therefore, for each m one obtains the untrue identity 1 
1 = -1. 



-2m+l 
2m+l ' 



which in the limit of m 



(20) 

oo becomes 



For the even number of quTVits and the multiports N — 2m (m > 2) we choose the following phases: 



2m-l ' • 



2m- 



Y-) = 0, 0jv = (0, 0, 0) = 0'. One easily finds that 



Eqm(< 
1 
2m 



1 0') = E QM (4>',- 



= ... = E, 



QM\ 



,0',' 



(2m — 1) exp 



,2m- 1 
' 2m - 1 ' 



exp i 



(2m - l) 2 
2m - 1 



(21) 



and 



Eqm(4>, --,0) 
1 



2m 
1 

2m 



(2m — 1) exp 
(2m — 1) exp 



2m 



2m- 1 
2m 

-i 7T 

2m- 1 



.2m(2m+ 1) 

exp | l 7T 

1 5 2m + 1 



Applying the same reasoning as above one has 1 = ^ (2m — 1) exp ^~» 2 m"i ' 
one has a contradiction, which in the limit of m — > oo becomes 1 = — 1. 



(22) 

Again, for each m > 2 



V. CONCLUSIONS 



We have derived the series of paradoxes in which the number of observers and the dimension of the Hilbert 
space describing each subsystem is the same. The derivation of the paradoxes relies on the perfect correlations 
observed in the system but the resulting contradiction between quantum mechanics and local realism manifests 
itself in the fact that local realism predicts a certain perfect correlation whereas quantum mechanics does not. 

An interesting feature of the paradoxes is that they naturally split into two parts: the even and the odd 
number of observers (input ports). In each part the final contradiction has different numerical values. All 
paradoxes do not vanish with the growing dimension N of the Hilbert space of each quiVit. 

Furthermore, for the case of two observers one cannot derive the GHZ paradoxes with the method presented 
here (for this case we have so called Hardy paradox but for non-maximally entangled state p5[). 
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In conclusion we state that the multiport beam splitters, and the idea of value assignment based Bell numbers, 
lead to a strikingly straightforward generalisation of the GHZ paradox for N entangled quiVits. These proper- 
ties may possibly find an application in future quantum information and communication schemes (especially as 
GHZ states are now observable in the lab 
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